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Abstract. We study the Green function gr r for the Laplace operator on the quotient of 
the hyperbolic plane by a cofinite Fuchsian group T. We use a limiting procedure, starting 
from the resolvent kernel, and lattice point estimates for the action of T on the hyperbolic 
plane to prove an "approximate spectral representation" for gr r . Combining this with 
bounds on Maafi forms and Eisenstein series for T, we prove explicit bounds on gr r . 

1. Introduction and statement of results 

(N ■ 

The hyperbolic plane H is the unique two-dimensional, complete, connected and simply connected 
Ricmannian manifold with constant Gaussian curvature —1. We identify H with the complex 
upper half-plane; this gives H a complex structure. In terms of the standard coordinate z = x + iy, 
| the Ricmannian metric is 

, ^ , dz dz dx 2 + dy 2 



and the associated volume form is 



(3z) 2 y 2 

i dz A dz dx A dy 



, Instead of using the geodesic distance r(z,w) on H directly, we use the more convenient function 

u(z,w) = coshr(z,u>) 

_ 1 \z-w\ 2 

(NJ ' 2(3?z)(3?w)' 
> ' 

; Let A denote the Laplace-Beltrami operator on H, given by 

Si A=y'(d> + %). 

| The Green function for A is the unique smooth real-valued function gr H outside the diagonal 

£SJ ■ on H x H satisfying 

gr H (z, w) = — log \z — w\ + 0(1) asz4»), 

/\ • Agr H ( ,w) — S w for all w G H, 

gr H (z, w) = 0(u(z, w)~ x ) as u(z, w) —> oo, 

where A is taken with respect to the first variable. It is given by 

gr H (z,u>) = -L(u(z,w)), 

where 

L( u ) = -LlogH±l. (1.1) 

47T U — 1 
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The group SL 2 (R) acts on H by isometries. Under the identification of H with the complex 
upper half-plane, this action on H is the restriction of the action on P 1 (C) by Mobius transfor- 
mations. Elements of SL2(R) \ {±1} are classified according to their fixed points in P 1 (C) as 
elliptic (two conjugate fixed points in P 1 (C) \ P 1 (R)), parabolic (a unique fixed point in P 1 (R)), 
and hyperbolic (two distinct fixed points in P 1 (R)). This terminology also applies to conjugacy 
classes. 

A Fuchsian group is a discrete subgroup of SL 2 (R). A cofinite Fuchsian group is a Fuchsian 
group T such that the volume of T\H with respect to the measure induced by /z H is finite. We 
will exclusively consider cofinite Fuchsian groups, and for such a group T we write 

voir = / A*h- 
Jr\H 

Remark. We define integration on T\H in a stack-like way so the above integral is l/#(rn{±l}) 
times the integral over T\H viewed as a Riemann surface. This implies that if / is a T-invariant 
function on H and T' is a subgroup of finite index in T, then 

/ / Mh = (r : r') • / /Mh- 
Jr'\H Jr\H 

Furthermore, this definition justifies the method of "unfolding" : if / is a smooth function with 
compact support on H and F is the function on T\H defined by 

then 

/ *V*H = / /MH- 

Jr\H Jh 

Let T be a cofinite Fuchsian group. The restriction of the Laplace operator A to the space 
of smooth and bounded T-invariant functions on H can be extended to an (unbounded, densely 
defined) self-adjoint operator on the Hilbert space L 2 (r\H), which we denote by Ar- 

The operator Ar is invertible on the orthogonal complement of the constant functions in the 
following sense: there exists a unique bounded self-adjoint operator Gr on L 2 (r\H, /x H ) such that 
for all smooth and bounded functions / on T\H the function Grf satisfies 

A r G r / = f - ^- [ / Mh and f Gt/Mh = 0. 
vol r Jr\n Jr\n 

There exists a unique function gr r on H x H that is T-invariant in both variables separately, 
satisfies gr r (z,u>) = gr r (w,z), is smooth except for logarithmic singularities at points of the form 
(2,72), and has the property that if / is a smooth and bounded r-invariant function on H, then 
the function Grf is given by 

Grf(z) = / gr r {z,w)f(w)n H (w). 
Jwer\u 

The function gr r is called the Green function of the Fuchsian group T. 

In this paper, we study gr r quantitatively, with the goal of obtaining explicit upper and lower 
bounds. One result that can be stated without introducing too much notation is the following. 

Theorem 1.1 (corollary of Theorem 4.1). Let r be a cofinite Fuchsian group, let Y be a compact 
subset of To\H, and let 5 > 1 and n > be real numbers. There exist real numbers A and B 
such that the following holds. Let T be a subgroup of finite index in Tq such that all non-zero 
eigenvalues of — Ar are at feast n. Then for all z, w e H whose images in r \H lie in Y , we have 

A < gr r (z, w)+ ( L ( M ( Z ' - L (' 5 )) ^ B - 

u(z,~fw)<5 
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Let us give a very brief overview of the article. In Section 2, we collect known results about 
Fuchsian groups. Most importantly, we make use of the hyperbolic lattice point problem and 
the techniques used to attack this problem. In Section 3, we use these results, together with 
a construction of the Green function involving the resolvent kernel, to "sandwich" the Green 
function gr r (with the logarithmic singularity removed) between two functions that, unlike gr r 
itself, admit spectral representations. In Section 4, we bound these functions in a way that lends 
itself to explicit evaluation. As an example, we find explicit constants A and B as in Theorem 1.1 
in the case where T = SL 2 (Z), T C T is a congruence subgroup, 5 — 2 and Fo is the compact 
subset of To\H corresponding to the points z in the standard fundamental domain of SL2(Z) such 
that 3z < 2. In Section 5, we use the bounds given by Theorem 1.1 to deduce bounds on gr r (z, w) 
in the case where Y is obtained by cutting out discs around the cusps of Y and where z, w or both 
are in such a disc. Finally, a number of bounds on Legendre functions that we will need have been 
collected in an appendix. 

In a forthcoming paper, we will use the results in this article to obtain explicit bounds on the 
canonical Green function of a modular curve X. This function is defined similarly to gr H , using 
the canonical (1,1) -form on X instead of fi H . It plays a fundamental role in Arakclov theory; 
see Arakelov [1] and Faltings [6]. Explicit bounds on canonical Green functions of modular curves 
are relevant to the work of Edixhovcn, Couveignes et al. [3] and the author [2] on computing 
two-dimensional representations of the absolute Galois group of Q that are associated to Hecke 
eigenforms over finite fields. 

This article may be compared with earlier work of Jorgenson and Kramer on bounding canon- 
ical and hyperbolic Green functions of compact Riemann surfaces [11, especially Theorem 4.5]. 
Jorgenson and Kramer consider compact Riemann surfaces X of genus at least 2, which can be 
obtained as X = Y\H for a cofinitc Fuchsian group Y without elliptic and parabolic elements. 
They obtain bounds on the hyperbolic Green function by comparing it to the heat kernel on X. 
Our method, too, starts with comparing the Green function with a kernel that can be obtained as 
a sum over elements of T, but the subsequent arguments are rather different. Let us note some of 
the differences. First, we allow arbitrary cofinite Fuchsian groups, which is the natural setting for 
modular curves. Second, the procedure that we apply in §3.1 to construct the Green function as 
a limit of a family of kernels K a for a — > 1 leads to bounds that are independent of the specific 
family. We take K a to be the resolvent kernel with parameter a —¥ 1, but the heat kernel with 
parameter t — > oo could have been used with the same result; sec [2, §11.5.2]. Finally, our bounds 
are much easier to make explicit than those in [11]; this is illustrated in §4.3. 

Acknowledgements . Part of this paper was written during a stay at the Max-Planck- Institut fiir 
Mathematik in Bonn; I am grateful for its hospitality. I thank Ariyan Javanpeykar for comments 
on an earlier version. The computations outlined in §4.3 were carried out using PARI/GP [14]. 

2. Tools 

2. 1 . Cusps 

Let T be a cofinite Fuchsian group. The cusps of Y correspond to the conjugacy classes of non- 
trivial maximal parabolic subgroups in Y. Every such subgroup has a unique fixed point in P (R). 
For every cusp c wc choose a representative of the corresponding conjugacy class and denote it 
by T c . 

Let c be a cusp of Y. We fix an element er c e SL 2 (R) such that er c oo <G P X (R) is the unique 
fixed point of T c in P 1 (R) and such that 

{±iK- 1 r c a c = {±i}{(^) |&ez}. 

Such a <7 C exists and is unique up to multiplication from the right by a matrix of the form ± ( J ^ ) 
with x € R; see Iwaniec [10, § 2.2]. We define 

q c :U^ C 

z i ^ exp(27riCT^ 1 z) 



3 



and 

y c :H^ (0,oo) 



k>g\q c (z)\ 



z i V xscr c z = — 

ZTT 

For all 7 e T, we write 



C c ( 7 ) = \c\ ifa"Vc 

Then we have 



a b 

c d 



r c = { 7 er|c c ( 7 ) = o}. 

It is known that the set {C c (j) | 7 e T, 7 ^ T c } is bounded from below by a positive number, and 
that if e is a real number satisfying the inequality 

< e < min C c ( 7 ), (2.1) 

then for all z e H and 7 e T one has the implication 

y c (z) > 1/e and y c (7z) > 1/e 76 T c . 
For any e satisfying (2.1), the image of the strip 

{x + iy I < x < 1 and y > 1/e} C H 

under the map 

H ^> H > T\H 

is an open disc D c (e) around c, and the map q c induces a chart on T\H identifying D c (e) with 
the punctured disc {z € C | < \z\ < exp(— 2ir/e)}. A compactification of T\H can be obtained 
by adding a point for every cusp c in such a way that q c extends to a chart with image equal to 
the disc {z € C | \z\ < exp(— 2n/e)}. Let D c (e) denote the compactification of D c (e) obtained by 
adding the boundary dD c (e) in T\H and the cusp c. 

Remark. Let us fix a point w E H and write for the stabilisator of w in T. The behaviour of 

gr r (z,w) as z — > w is 

gr r (z, w) = log \z — w\ as z — > w. 

27T 

Furthermore, the behaviour of gr r (z, w) as z moves toward a cusp c of T is 

gr r (z,w) = -i-log2/ c (z) + 0(l) as y c (z) 00. 
2.^. T/ie Selberg-Harish- Chandra transform 

Let 3: [l,oo) — > R be a smooth function with compact support. The invariant integral operator 
attached to g is the operator T g defined on smooth functions /: H — >• C by 



(T g f)(z)= [ g(u(z,wj)f(w)n H (w). 
J wen 



1 wen 

The Laplace operator A commutes with all such operators T g ; sec Sclbcrg [16, pages 51-52] or 
Iwaniec [10, Theorem 1.9]. In fact, every eigenfunction of A is also an eigenfunction of all invariant 
integral operators, and conversely; see Selberg [16, page 55] or Iwaniec [10, Theorems 1.14 and 1.15]. 
The relation between the eigenvalues of A and those of T g is given by the Selberg-Harish- Chandra 
transform of g. This is a holomorphic function h defined by the following property. Let /: H — > C 
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be an eigenfunction of —A with eigenvalue A = s(l — s). Then / is also an eigenfunction of T g , 
and the eigenvalue depends only on A; we can therefore define h(s) uniquely such that 

-A/ - s(l - a) J => T g f = h(s)f. (2.2) 

In particular, taking / = 1, we see that 

/oo 

g(u)du. (2.3) 

The Selberg-Harish-Chandra transform can be identified with the classical Mehler-Fock trans- 
form, defined as follows (see Iwaniec [10, equation 1.62']): 

/oo 
g{u)P s -i(u)du. (2.4) 

Here P v is the Legendre function of the first kind of degree v; see Iwaniec [10, equation 1.43] or 
any book on special functions, such as Erdelyi et al. [4, §3.6.1]. The function g can be recovered 
from h, and we call g the inverse Selberg-Harish-Chandra transform of h. 

The identity (2.2) holds more generally than just for smooth functions g with compact support; 
see Selberg [16, pages 60-61]. It will be enough for us to state a slightly weaker, but more convenient 
sufficient condition (cf. Selberg [16, page 72] or Iwaniec [10, equation 1.63]). Let e > and /? > 1, 
and let h be a holomorphic function on the strip {s G C | -e < Ss < sucn that h(s) = h(l — s) 

and such that s i->- |/i(s)||s(1 — s)| /3 is bounded on this strip. Then the inverse Selberg-Harish- 
Chandra transform g of h exists, and (2.2) holds for the pair (g, h). 

2. 3. Spectral theory of the Laplace operator for Fuchsian groups 

Let L be a cofinite Fuchsian group. The spectrum of — Ar on L 2 (F\H) consists of a discrete part 
and a continuous part. 

The discrete spectrum consists of eigenvalues of — Ar and is of the form {Xj}°^ with 

= Ao < Ai < A2 < . . . , Xj — > 00 as j — > 00. 

Let {cfij}jjL be a corresponding set of eigenfunctions; these are called automorphic forms of Maafi 
( of weight 0). We may and do assume that they are orthonormal with respect to the inner product 
on L 2 (F\H). For each j > 0, we define s <G C by 

^3 = s i(l - s j)j 

with Sj G [1/2, 1] if Xj < 1/4. For Xj > 1/4, the Sj are only determined up to Sj -H- 1 — Sj. 

The continuous part of the spectrum of — Ar is the interval [1/4, 00) with multiplicity equal 
to the number of cusps of T. In particular, the continuous spectrum is absent if T has no cusps. 
The continuous spectrum does not consist of eigenvalues, but corresponds to "wave packets" that 
can be constructed from non-holomorphic Eisenstein series or Eisenstein-Maafi series, introduced 
by Maafi in [13]. These series are defined as follows: for every cusp c of Y the series 

E c (z, s) = (v^ z )Y (z E H, s e C with 3£s > 1) 

7er c \r 

converges uniformly on sets of the form K x {s £ C | Uis > 5} with K a compact subset of H and 
5 > 1. In particular, E c (z, s) is a holomorphic function of s. 

A crucial ingredient in the spectral theory of automorphic forms is the meromorphic contin- 
uation of Eisenstein series, due to Selberg [17]. For proofs of this meromorphic continuation and 
of the other properties of the Eisenstein series that we use, wc refer to Hcjhal [8, Chapter VI, 
§11]. Different constructions of the meromorphic continuation can be found in Faddeev [5, §4], 
Hejhal [8, Appendix F] or Iwaniec [10, Chapter 6]. 

The meromorphic continuation of the Eisenstein series takes the following form. The functions 
E c (z,s) can be continued to functions of the form E c (z,s) — H(z,s)/G(s), where H is a smooth 
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function on T\H x C and both G and H are entire functions of s. These meromorphic continuations 
have a finite number of simple poles on the segment (1/2, 1] and no other poles in {s £ C | 3?s > 
1/2}, and satisfy a functional equation, which we will not write down. For all s £ C that is not a 
pole, the function z H> E c (z, s) satisfies the differential equation 

-A r E c ( ,s) = s(l-s)E c ( ,s). 

For s £ C with 3ts = 1/2, the Eisenstein-Maafi series E c ( , s) are intcgrablc, but not square- 
integrable, as functions on T\H. In contrast, the "wave packets" mentioned above are square- 
integrable. They figure in Theorem 2.1 below, which is a fundamental result in the theory of 
automorphic forms. 

In the following theorem, and in the rest of the article, we will consider integrals over the 
line JJs = 1/2. For this we need an orientation on this line; we fix one by requiring that the map 
1 1-> 1/2 + it from R with the usual orientation to the line 9?s = 1/2 preserves orientations. 

Theorem 2.1 (see Iwaniec [10, Theorems 4.7 and 7.3]; cf. Faddeev [5, Theorem 4.1]). Every 
smooth and bounded T-invariant function f : H — > C has the spectral representation 

oo . 

/(*) = !>&(*) + E T7 / b c (a)E c (z,a)ds, (2.5) 

j=0 c 4 ™ J*t'=l/2 

where c runs over the cusps of T and the coefficients bj and b c (s) are given by 

bj= / Mmh and b c (s)= / fE c ( ,s)m h - 
Jr\n Jr\n 

The right-hand side of (2.5) converges to / in the Hilbert space L 2 (r\H). If in addition the smooth 
T-invariant function A/: H — > C is bounded, the convergence is uniform on compact subsets of H. 

With regard to the spectral representations provided by this theorem, the effect of the oper- 
ator Gr from the introduction is as follows: if / has the spectral representation (2.5), then Grf 
has the corresponding spectral representation 

*/W - " | £*M - Y. 55 L 1/2 ■)*■ <«» 

(Note the absence of the eigenvalue A = 0.) 

There is an analogous result (Theorem 2.2 below) for functions on H x H that arc of the 
form X) 7 er 9( u ( z > 7 w ))j where g: [l,oo) — > R is the inverse Selberg-Harish- Chandra transform of 
a function h as at the end of § 2.2. The result involves a type of convergence that we now explain. 
Let A be a filtered set, and let {K a } a£ A be a family of continuous functions on T\H x T\H that 
are square-integrable in the second variable. If K is a function such that for all compact subsets C 
of F\H we have 

liml sup \K a (z, w) - K(z, w)\ + sup / \K a (z, w) — K(z, w)| 2 /i H (w) ) = 0, 

a ^ A \z,wec zecJ w er\H J 

we say that the family of functions {K a } ae A converges to K in the (Lj^ c ,L 2 n LJ^ C ) -topology . In 
other words, this condition means that the family converges uniformly on compact subsets of Hx H, 
and also with respect to the L 2 -norm in the variable w, uniformly for z in compact subsets of T\H. 

Theorem 2.2 (see Iwaniec [10, Theorem 7.4]). Let g: [l,oo) — >• R be the inverse Selberg-Harish- 
Chandra transform of a function h as at the end of § 2.2. Then the function 

K g : H x H > R 

(z,w) i — ► ^2g(u(z,^w)) 

is T-invariant with respect to both variables and admits the spectral representation 
oo _ If 

K g (z,w) = y2h(s j )cf> j (z)(f> j (w) + y2 T ^ h( S )E c (z,s)E c (w,s)ds, (2.7) 

where the expression on the right-hand side converges to K g in the (Lg^L 2 n Lj^ c ) -topology. 
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2.4- A point counting function 

We fix a real number U > 1, and define 

9u- [l,oo) — > R 

fl ifw<[/; ( 2 - 8 ) 
M lO if u > 17. 

From (2.4) and the formula for J* P u (w)dw found in Erdelyi et al. [4, §3.6.1, equation 8], we sec 
that the Selberg-Harish-Chandra transform of gu is 



h v {s) = 2tt^U 2 -IP~\(U). (2.9) 

Here Pjf is the Legendre function of the first kind of degree v and order /1; see [4, § 3.2]. 

Now let r be a cofinite Fuchsian group. We introduce the following point counting function. 
For any two points z, w in H and any U > 1, we denote by Nr(z, w, U) the number of translates 
of w by elements of T lying in a disc around z of radius r given by cosh(r) = U, i.e. 

N r (z, w, U) = #{7 e T I u(z, 7w) < U} 

= Y / 9u(u(z,w)). ( 2 - 10 ) 

This is T-invariant in z and w separately. 

Lemma 2.3. Let U e [1, 3], and Jet s <= C be such that s(l - s)(t/ - 1) e [0, 1/2]. Then /it/(s) is 
a rcai number satisfying 

(4tt - 8)(I7 - 1) < hu(s) < 8(U - 1). 
Proof. This follows from (2.9) and Lemma A.l. □ 
2.5. Bounds on eigenf unctions 

The convergence of the spectral representation (2.7) can be deduced from suitable bounds on the 
function 

$ r :H x [0,oo) — ► [0,oo) 

(z,X)^ E l^)| 2 + £^-L=i/2 \Ec(z,a)\ a d.. ( 2J1 ) 

j-.XjKX c J a (l-s)<A 

We will prove a bound on $r which holds uniformly for all subgroups T of finite index in a given 
Fuchsian group r . This will give a similar uniformity in Section 4. 

Lemma 2.4. Let T be a cofinite Fuchsian group. Then the function <J>r(z, A) satishes 

$r(z, A) < ^ JV r (;e,;8,17)A for aii z e H and aii A > 1/4. 
Proof. Let z e H and A > 1/4. We put 

U = i + ±e(i,3). 

From Bessel's inequality one can deduce (see Iwaniec [10, §7.2]) that 

E IM^Wf + Er? L s =i/ 2 |M*)£c(*,*)| 2 <k< / iVr( Z ,^C/) 2 MH H- 

J: A,<A c 47r *- / .(i_ 8 /<A J - er \ H 

Using the definition (2.11) of $r and the bound hu(s) > (2n — 4)/A given by Lemma 2.3, we 
deduce 

*r(«,A) < A / AHz^LOVhH- 
(2tt - 4)^ y wer \H 
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We rewrite the integral on the right-hand side by partial "unfolding" as follows (cf. Iwaniec [10, 
page 109]): 



/ N r (z,w,U) 2 fj, H (w)= V / g u {z,-f'w)g u {-fz,-f'w)^ u {w) 
Jwer\n 77 , gr Ju.er\H 

= gu{z,w)gu(jz,w)n H (w). 
7er J wen 

The last integral can be interpreted as the area of the intersection of the discs of radius r around 
the points z and jz of H, where cosh r = U . By the triangle inequality for the hyperbolic distance, 
this intersection is empty unless 

u(z,jz) < cosh(2r) = 2U 2 - 1; 
furthermore, the area of this intersection is at most 2n(U — 1) = n/X. From this we deduce that 

N T {z,w,U) 2 ^{w) < ^N r (z,z,2U 2 -l) 



I 

Ju 



™er\H ^ 

Since 2U 2 — 1 < 17, this proves the lemma. □ 
2.6. The hyperbolic lattice point problem 

Let T be a cofinite Fuchsian group. The hyperbolic lattice point problem for T is the following 
question: what is the asymptotic behaviour of the point counting function Nr(z, w, U) from (2.10) 
as U — > oo? In contrast to the Euclidean analogue of this question, about the number of points 
in Z 2 lying inside a given circle in R 2 , no elementary method is known to even give the dominant 
term. The difficulty is that for circles in the hyperbolic plane of radius tending to infinity, the 
circumference grows as fast as the enclosed area: the circumference of a circle of radius r equals 
27rsinh(r), and the area of a disc of radius r equals 27r(cosh(r) — 1). In spite of this difficulty, good 
estimates for Np(z, w, U) can still be found, namely using spectral theory on T\H. 
The strategy is to take suitable functions 

9u,9u'- 1 1 ' 00 ) -> R 

with compact support, and to define functions Ky and KjJ on H x H, invariant with respect to 
the action of T on each of the two variables, by 

K±(z,w) =^g^{u{z,-iw)). 

This sum is finite because the functions gfj have compact support. We take the functions g^j such 
that for all z, w e H and U > 1, we have the inequality 

K v {z,w) < N r (z,w,U) < K+{z,w). (2.12) 

Provided the Selberg-Harish-Chandra transforms h v of g v satisfy the conditions of Theorem 2.2, 
the functions Ky have spectral representations 

oo _ . 

^(^) = E^(4W+E« /_ h±(s)E c (z,s)E c (w,s)ds. (2.13) 

j=0 c JSRs— 1/2 

These spectral representations can then be used to find the asymptotic behaviour of Nr(z,w, U) 
as U — >• oo. 
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A reasonable choice at first sight would be to take for both g^ and the function gu defined 
by (2.8), so that the inequalities in (2.12) become equalities. Unfortunately, the Selberg-Harish- 
Chandra transform hjj of gu does not decay quickly enough as \$ss\ — > oo to give a spectral 
representation of Nr(z,w,U) as in Theorem 2.2. We take instead 

r 1 if 1 < u < U, 
9u(u) = \ ffi if (7< U <U, 
I if V < u 

and 

r 1 if 1 < u < T, 
9u(u) = \ if T <u<U, 

I if U < u 

for certain T, V, depending on [/, with 1 < T < U <V; see Iwaniec [10, Chapter 12]. Using (2.4), 
we obtain 

U-u 



V-U 



du, 



P s _ 1 (u)y—jjdu-2ir J P a _i(«) 

Ps-l(u)jJ-^du-2TT J P^^-^—du. 

Integrating by parts and applying the integral relation between the Legendre functions P v and P~ 2 
given in Erdelyi ct al. [4, §3.6.1, equation 8], we get 

, 9 {v 2 -i)p;\(v)-{u 2 -i)p;\{u) 

h U s ) = 27r t7 77 > 

(2 14) 

h - ( , 9 (c/ 2 -i)P s -i 2 i(^)-(r 2 -i)n- 2 i(r) 

M s ) = 2?r ^"7^ • 

The dominant term in (2.13) as U — > oo comes from the eigenvalue Ao = 0, corresponding to 
s = 1. It follows from (2.3) or the formula Pq 2 (u) = (u — l)/{2u + 2) that 

h+(l) = 2n(U-l)+n(V-U) and h u (l) = 2ir(U - 1) - n(U - T). (2.15) 

Let a real number 8 > 1 be given. We fix parameters a + , a~ , /?+ and /3~ satisfying 

a±e (0,1/2), /3±>0, /T<— -. (2.16) 

o + l 

We choose T and V as functions of U as follows: 

T{U) = U-P~ U- x - a ~(U 2 - 1), U(C7) = J7 + /3+?7- 1 - a+ (C/ 2 - 1). (2.17) 

The last inequality in (2.16) ensures that iiU>8, then T(U) > 1. 

For later use, we will keep the parameters and ^ variable for greater flexibility. To obtain 
the best known error bound in the hyperbolic lattice point problem, the right choice is a* = 1/3, 
so that 

V -U ~ p+U 2/3 and U -T ~ fi-U 2/3 as U -> oo. 
This choice leads to the following theorem. 

Theorem 2.5 (Hubcr [9, Satz B], Patterson [15, Theorem 2], Sclbcrg; see Iwaniec [10, Theo- 
rem 12.1]). Let r be a cofinite Fuchsian group. For all z,w <G H, the point counting function Nr 
satisfies 

N r (z,w,U)= 2 ''^ S + I! hWfcW' + 0{U 2 ' 3 ) asC/^oo, 

with an implied constant depending on T and the points z and w. 
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In particular, since \4>o\ 2 is the constant function 1/volr, this shows that 

m 27r(f/-l) 

N r (z,w,U) ~ — ^— ; as£/^oo. 

voir 

Since 27r([7 — 1) is the area of a disc of radius r with coshr = U, Theorem 2.5 implies that this 
area is asymptotically equivalent to the number of lattice points inside the disc times the area of 
a fundamental domain for the action of T, which is the intuitively expected result. 

3. An approximate spectral representation of the Green function 

Let r be a cofmite Fuchsian group. The Green function of T formally has the spectral representation 

? 00 i i r i 

gr r (z,w) =-J^T M z )M w )~^2t^ / 771 ~:E c (z,s)E c (w,s)ds. 

The problem is that this expansion does not converge. Neither should one be tempted to write the 
Green function by "averaging" gr H as a (likewise divergent) sum 

gr r (z,u;) = ^gr H (z, 7 w). 

However, both of these divergent expressions have at least some value as guiding ideas for what 
follows. In fact, we will bound gr r (z, w) by means of certain functions Rp s (z,w), defined in (3.5) 
below, that reflect the above formal spectral representation of gr r . 

3.1. A construction of the Green function using the resolvent kernel 

We will give a construction of the Green function of T using the family of auxiliary functions 

g a : (l,oo) — > [0,oo) 

U I > ^Q a -l{u) 

for a > 1, where Q v is the Legendre function of the second kind of degree v\ see Erdelyi et al. [4, 
§3.6.1]. By [4, §3.6.2, equation 20], we have 

Qo(u) = -log— ^ 

which shows that g\ equals the function L from (1.1). By (2.4) and [4, §3.12, equation 4], the 
Selberg-Harish-Chandra transform of g a is 

/oo 
P s -i{u)Q a -i(u)du 

1 



(a — s)(a — 1 + s) 
1 



s(l - s) +a(a- 1)' 
Given a real number a < 1/2, we consider the strip 

S a = {s e C | a < 5fts < 1 - a}. (3.1) 
Lemma 3.1. (a) For all a, b > 1 and all a e (1 — min{a, b}, 1/2), the function 

I II 1 2 

s i y \h a (s) - hb(s)\\s(l - s)\ 

is bounded on 5 CT . 

(b) Let a e (0, 1/2). There exist real numbers {C a fi,a)a,b>i, with C a ,b,a — > as both a and b tend 
to 1, sucii that 

I I I I 2 

\h a {s) - h b (s)\ < CaM,<r\s(l - s)\ ~ forallseS a . 
Proof . Both claims are easily deduced from the expression 

6(6- l)-o(o- 1) 



K{s) - h b (s) = 



(a - s)(a - 1 + s)(b - s)(6+ 1 - s) 
Details are left to the reader. □ 
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For alia > 1, the sum X^ 7 gr 9a{u(z, jw)) converges uniformly on compact subsets of H x H not 
containing any points of the form (z, jz) and defines a continuous function that is square-integrable 
in each variable; see Fay [7, Theorem 1.5]. We can therefore define 

K T a .{{z,w) eHxH|z^r«)HR 

0, w) i— > 9a(u(z, t«;)) - Co, ( 3 - 2 ) 

where 

27T f°° ( ^ 1 i, m 1 

Ca = — j- / 9a[u)du = — — ft a (l) = — ■ — . 

voir Ji voir voir a (a — 1) 

The constant c a is such that the integral of K% over T\H with respect to each of the variables 
vanishes. Up to this constant, is the resolvent kernel with parameter a. 

It is known that the resolvent kernel admits a meromorphic continuation in the variable a. 
The following result can be interpreted as the statement that — gr r is the constant term in the 
Laurent expansion of the resolvent kernel at a = 1. 

Proposition 3.2. The family of functions {— K^} a>1 converges to the Green function gr r in the 
(L-,L 2 nL-)-topoJog 7 . 

Proof. It follows from Lemma 3.1(a) that for all a,b > 1, the function g a ~g b satisfies the conditions 
of Theorem 2.2. The function 

( K a - K b )(z,w) = ^2(g a (u(z,jw)) - g b (u(z,jw))) -c a + c b 
therefore has the spectral representation 

oo 

(Kl K b r )(z,w) = J2(h a ( Sj ) - h b ( S3 ))<f> 3 (z)^(w) 

3=1 (3.3) 
If 

+ V— / (h a (s) - h b (s))E c (z,s)E c (w,s)ds, 
4tti J Us=1/2 

where the right-hand side converges to — in the (L^ C ,L 2 n L^ c )-topology. (Note that the 
eigenvalue Ao = has disappeared because of the definition of c a .) In particular, — extends 
to a continuous function on H x H that is T-invariant with respect to both variables. 

We claim that {K^ — if^} a ,fc>i converges to in the (Lg^, L 2 n L^ c )-topology as a, b \ 1. In 
particular, this implies that {K^ } a> i converges to a symmetric continuous function on T\H x T\H 
that is square-integrable with respect to each variable separately. We fix ct G (0, 1/2) be such that 
the spectrum of — Ar is contained in {0} U [a(l — a), oo). 

First we show that {K^ — if^} a ,&>i converges to zero uniformly on compact subsets of H x H. 
Lemma 3.1(b) implies 

oo 

K - K T b \{z,w) < J2 K( Sj ) h b ( Sj )\ ■ |^(z)^H| 

3=1 

+ V— / \h a {s)-h b (s)\-\E c {z,s)E c {w,s)\ds 

(oo 
£( Sj (i- Sj 0r 2 |<M^MI 
3 = 1 

+ Y,J- [ (s(l-s))~ 2 \E c (z,s)E c (w,s)\ds 

By the Cauchy-Schwarz inequality and Lemma 2.4, the right-hand side converges to uniformly 
on compact subsets of H x H, as claimed. 
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Next we show that {K^ — -K^} a ,i»i converges to zero with respect to the L 2 -norm in the 
variable w, uniformly for z in compact subsets of H. From (3.3), Plancherel's theorem and 
Lemma 3.1(b), we deduce 



/ 



«>er\H 



K a - Kb\ 2 {z,w)ii H (w) = ^2\h a (sj) - h b (sj)\ \<t>j{z)[' 



+ y~\-T—- f \K{s) - h b (s)\ 2 \E c (z,s)\ 2 ds 



— ™» Jffts=l/2 
/ oo 

< C a,b,a - S j ))~ 4 \<p j (z)\' 1 



+ Ep / (s(l-s))- 4 \E c (z,s)\ 2 ds) 



for real numbers C 0j b i(T with C a .b,a — > as a, b \ 1. By Lemma 2.4, the other factor on the right- 
hand side is bounded on compact subsets of H. This implies that the right-hand side converges 
to uniformly on compact subsets of H as a, b \ 1, as claimed. 

The defining property (2.2) of the Selberg-Harish-Chandra transform implies that if / is a 
smooth, bounded, L-invariant function on H, with spectral representation (2.5), then 

/ Kl(z,w)f(w)fj, Ii (w) = SZbjhaisjWjiz) + I b c (s)h a {s)E,(z 7 s)ds. 

Jwer\H ~[ ( 47r * Js? s =i/2 

Taking the limit, using the L 2 -convergence that we just proved and applying (2.6), we get 
lim Kl(z,w)f(w)n H (w) = lim / K%(z, w)f(w)n Ii {w) 

= -Grf(z) 



J 

Jwt 



gr r (z,w)f(w)(j, a (w). 

lwer\n 

Since the set of smooth and bounded functions is dense in L 2 (r\H), this proves that the limit of 
the convergent family of functions {K^ } a >i equals — gr r . □ 

3.2. Proof of the approximate spectral representation 

We now exploit the estimates for the hyperbolic lattice point problem given in § 2.6. We choose 
parameters a ± and /3 ± satisfying (2.16). Using these, we define functions T(U), V(U), g^(u), 
hy(s) and K^(z,w) as in §2.6. Furthermore, we fix a real number S > 1. We write T as the 
disjoint union of subsets Ur t s(z, w) and Ar t s{z,w) defined by 

U r ,g(z,w) = {7 E T I u(z,jw) < 6}, 
Ar,s{z,w) = {7 e L I m(z,7w) > 5}. 



We define 



1 r°° h ± (i\ 

^ {S) =2^J S TP^[ dU f ° r < ^ s < ( 3 - 4 ) 
R± 5 (z,w) =J2 I s t (sj)<t>j(z)hM + Y,T1 / lt(s)E t (z,8)E c (w,s)ds, (3.5) 

+ 1 ( P + 1 6+1 \ - 1 ( ^ j.1 S+1 \ f*n 

q ^ = ^{2c7^- l0g —)> q ^ = -^{2^ +l0g —)- (3 ' 6) 
The intuition behind the following theorem is that although the Green function gr r does not 
admit a spectral representation, it can be bounded (after removing the logarithmic singularity) 
by functions that do admit spectral representations. The terms qp s below correspond to the 
eigenvalue 0, while the terms s (z, w) correspond to the non-zero part of the spectrum. 
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Theorem 3.3. Let T be a cofinite Fuchsian group. For all S > 1 and for every choice of the 
parameters and ^ satisfying (2.16), the Green function of T satisfies the inequalities 

-q+ s - R+ S (z,w) <gr r (z,w)+ Y i L ( u ( z ^w)) - L(5)) < -q~ s - R~ s {z,w). 

76n r ,5(2,i«) 

Proof. For any U > 5, the inequality (2.12) implies that the number of elements 7 e Ar,s(z,w) 
with u(z 7 jw) < U can be bounded as 

A(U) < #{7 e A r , 4 (*, w) I u(z, jw) <U}< B(U), (3.7) 

where the functions A, B: [<5, 00) — > R are defined by 

A(U)=K-(z,w)-#Il r , 6 (z,w) and B(U) = K+(z, w) - #n r ,«(*, w). 

The functions A and B are continuous and increasing. The estimates from §2.6 imply that they 
are bounded linearly in U as U — > 00, with an implied constant depending on the group T, the 
points z and u; and the functions T and V. 

Let {fi a }a>i, {5a}a>i and {K^} a> i be as in §3.1. For all a > 1, applying partial summation 
and (3.7) gives 

pOO _ />oo 

-/ g' a (U)A(U)dU < ]T 9a{u{z, 1W ))<- g' a (U)B(U)dU. 
Using the definition (3.2) of K^, we deduce the upper bound 

pOO c\ pOO 

K T a {z,w)< Y, 9a{u{z, 1W ))- g' a (U)B(U)dU - — g a (u)du. 
jeu rAz , w) J* volr J* 

The definition of B implies 

r-oo /*00 f-OO 

/ g' a (U)B(U)dU = g' a (U)K+(z,w)dU-#U rtS (z,w) g' a (U)dU 
Js Js Js 

f°° / 9tT \ 9-7T f°° 

= J e ^U)(K { z, w )-— { U-l)) d U + —J s g' a (U)(U-l)dU 
+ #n r , s (z,w)g a (5). 

Using integration by parts, we rewrite the second integral in the last expression as follows: 

/•oo poo />5 

/ g' a (U)(U-l)dU= g' a (U)(U-l)dU- g' a (U)(U - l)dU 

JS Jl Jl 

/oo pS 
9a (U)dU- J g' a (U)(U - l)dU. 

We can now rewrite our upper bound for K^(z,w) as 

K T a (z,w)< Y {9a(u(z,w)) - g a (6)) - jf g' a (U)[K+(z, w) - — (U - l))dU 

j\' a (U)(U - l)dU. 



7Gn r ,5(z,iii) 

+ 



2n '* 



voir 

Lemma A. 2 implies 
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^ / o \ a ~ 1 



and equality holds for a = 1. By the dominated convergence theorem, we may take the limit a \ 1 
inside the integrals. Together with Proposition 3.2, this leads to 

i poo O \ All 

gr r (z, w )+ J2 {L{u{z, 1W ))-L{5))>--j (K+(z,w) - —(U - 1)) 

1 (5+1 

+ ^: Iog — ■ 

In the integral, we insert the spectral representation (2.13) of if^, the formula (2.15) for hy(l) 
and the fact that \4>o\ 2 = 1/volr- We then interchange the resulting sums and integrals with the 
integral over U; this is permitted because the double sums and integrals converge absolutely, as 
one deduces from Lemma 2.4 and Theorem 2.5. This yields 

Finally, we note that 

roc -rr tj roo 

I —i du ^ + L u ^ a+du 

/3 + 



a+S a+ 

This proves the lower bound of the theorem. The proof of the upper bound is similar. □ 

Remark. The only inequality responsible for the fact that the inequalities in Theorem 3.3 are not 
equalities is (3.7). 

4. Explicit bounds 

4-.1. Bounds on hy(s) and I^(s) 

We keep the notation of § 3.2. In addition, we choose real numbers cr* such that 

<a + <cr+ < 1/2 and < a~ < a~ < 1/2. 

Let s be in the strip S a + defined by (3.1), and let p a +(u) be the elementary function defined 
by (A. 6) below. From (2.14), Corollary A. 6 and (2.17), we obtain 

WM |<2,g^^--' V -»l + < t/2 - 1 »l P ---' t/ »l 



< 2tt|s(1 - s)\ 
= 2tt|s(1 - s)\ 



v-u 

-5/4 Pa+ (V)+ Pa+ (U) 

v-u 

-5/4 {p a+ {V)+P«+(U))U 1+a + 

P+(U 2 -1) 



Similarly, for s € S a - , 

-5/4 {Pa-(U)+Pa-(T))U 1+C 



\hu(s)\ < 2n\s{l - s)\ 



p-{u 2 -i) 

Substituting this in the definition (3.4) of /, we obtain 

|/ 5 + ( S )|<^ + | S (l- S )|- 5/4 and |/,-(s)|<^-| S (l- S )|- 5/4 , (4.1) 



where 

n+ _ i [°° (pMv)+pMu))u 1+a+ , TT 



i r 



{p a -{U)+PaAT))U 1+a \ nT 

w^w a 



(4.2) 
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4-2. Bounds on gr r 

Theorem 4.1. Let F be a Fuchsian group. Let S > 1 and n G (0, 1/4] be reaJ numbers such that 
the spectrum of — Ar is contained in {0} U [77, 00). Let er+, er~, of 1 ", a - , /3~ be real numbers 
satisfying (2.16) and the inequalities 

0< a+ < er+ < 1/2, 0<a-<cr-<l/2 and a ± (l - ct ± ) < 77. 

Then the Green function of T satisfies the inequalities 

A<gr r (z,w) + E {L{u{z,jw))-L(5))<B for all z, w G H, 

7£ll ri {(z,ra) 

where 

A ~ qr ' s Us (2^4)2 \-^ +AV2 J 2 ' 

r „- n- E ^' 5/4 1 ,,/7T ^r(^^17) + jV r ( W , W ,17) 
S - + ^ (2^4)^ + 4V2 J 2 • 

Proof . In view of Theorem 3.3, we have to bound the absolute values of the functions i?p Az, w) 
from (3.5). Applying the triangle inequality and the Cauchy-Schwarz inequality, we see that 



where 5 + and 5* _ are defined by 

5±(.) = e i^-)ii^)i 2 + E 7- / i^wii^^s)! 2 ^. 

j=l c " , 3?s=l/2 

Let $r(^,A) be as in (2.11). Applying (4.1), we obtain (with $<I>r/<9A taken in a distributional 
sense) 

00 . . 

S ± {z)/Dt<Y. X T' i \^)\ 2 +llj- / (s(l-s)y 5/i \E c (z lS )\ 2 ds 



r i 00 5 Z" 00 

= rT 5/4 $rM/4) + A- 5 / 4 $r(*,A) + - A / A~ 9 / 4 <I>r(z, A)dA 

L JA=l/4 4J 1/4 

r /-CO 

= (tT 5/4 - 2 5 / 2 )$ r (^, 1/4) + 7 / A- 9 / 4 $r(^, A)dA. 

4 il/4 



The bound on $r(z, A) given by Lemma 2.4 implies 

SH*) < D±^L_tf r (*, z, 17) - 2^) . J + I X-^XdX 

= ^ ± (2^4F iVr ^^ 17) ( !Z X !+4 ^)- 
This proves the theorem. □ 

Proof of Theorem 1.1. Let the notation be as in the theorem; we may assume 77 < 1/4. We apply 
Theorem 4.1 to L, with parameters , and /J^ depending only on 77 and not on L. It is clear 
that the factor 1/volr occurring in the definition (3.6) is bounded by l/volr , and that Nr(z, z, 17) 
is bounded by Nr (z, z, 17). It remains to remark that Nr (z, z, 17) is bounded on Yq. □ 

The bounds given by Theorem 4.1 are easy to make explicit. First, the real numbers Df 
from (4.2) can be bounded in elementary ways or approximated by numerical integration. Second, 
a straightforward computation shows that for z = x + iy G H and 7 = ( a b d ) G SL2(R), we have 

l( f 2 | (b + (a — d)x — cx* 



u{z, 1Z ) = -y{a-cxY + ^ * ^ j + (cyY + (d + ex) 2 j . (4.3) 

This can be used for concrete groups F to find an upper bound on Nr(z, z,U) for U > 1. 
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4-3. Example: congruence subgroups o/SL 2 (Z) 

Let us consider the case where T = SL 2 (Z). We will make the bounds from Theorem 1.1 explicit 
for congruence subgroups T C SL 2 (Z). By our convention for integration on T\H if —1 £ T, we 
have 

1 1 7T 

< = -• 

voir vo1 S l 2(z) 6 

We choose 

6 = 2. 

Sclbcrg conjectured in [18] that the least non-zero eigenvalue Ai of — Ar is at least 1/4, and he 
proved that Ai > 3/16. The sharpest result known so far, due to Kim and Sarnak [12, Appendix 2], 
is that Ai > (25/64) (1 - 25/64) = 975/4096. We may therefore take 

r, = 975/4096. 

We now consider the point counting function N Sh2 ^(z 7 z, U) defined by (2.10) on a rectangle 
of the form 

R = {x + iy £ H | x min <x< x max , y min <y< y max } 

for given real numbers x m - m < x max and < y mm < J/max- The function z i-> -/VgL 2 (z)(<z, z, U) on R 
is clearly bounded from above by the number of matrices 7 = (" d ) £ SL 2 (Z) such that for some 
z £ R, the inequality 

u(z,"fz)<U (4.4) 

holds. We now show how to enumerate these matrices. We distinguish the cases c = and c / 0. 
We assume (after multiplying by —1 if necessary) that a = d = 1 in the first case, and that c > 
in the second case. The total number of matrices 7 as above is then twice the number produced 
by our enumeration. 

In the case c = 0, by (4.3), the inequality (4.4) reduces to 

l + \{b/yf<U. 

This implies 

l&l < 2WV2J7 - 2. 
In the case c > 0, it follows from (4.3) and (4.4) that 

|c| < V2~U/y min , 
-VZU + cx min <a< VZU + C.T max , 
-VZU - ca; max < d < VZU - cx min . 

Since c ^ 0, the cocfhcicnts a, c, d and the condition ad — be = 1 determine b. If 7 is a matrix 
obtained in this way, we compute the minimum of u(z, jz) for z £ R using (4.3) to decide whether 
there exists a point z £ R satisfying (4.4). 

Let Yq denote the compact subset in SL 2 (Z)\H which is the image of the rectangle 

{x + iy £ H I -1/2 < x < 1/2 and V3/2 <y<2}. 

This is the complement of a disc around the unique cusp of SL 2 (Z). Dividing this rectangle into 
100 x 100 small rectangles and bounding -/Vsl 2 (z) (^ z , U) on each of them as described above, we 
get 

N ro (z,z, 17) < 216 for all z £ Y a . 

Given this upper bound for iV SL2 ( Z ) (z, z, 17), some experimentation leads to the following values 
for the parameters: 

a+ = 0.0366, P + = 2.72, a+ = 0.306, 
a- = 2.96 • 10~ 3 , /3~ = 0.668, a' = 0.250. 

With these choices, a numerical calculation gives 

q+ s < 69.0, q^ 5 > -216, D+ < 18.5, DJ < 9.61. 

This implies the following explicit bounds on Green functions of congruence subgroups. 
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Theorem 4.2. Let T be a congruence subgroup of SL 2 (Z). Then for all z,w <G H whose images 
in SL 2 (Z)\H lie in Y , we have 



-2.87- 10 4 < gr r (z,w) + 



E 

u(z,'-yw)<2 



L(u(z,-/w)) - 



log 3 

4"7T 



< 1.51 • 10 4 . 



5. Extension to neighbourhoods of the cusps 

The bounds given by Theorem 4.1 do not have the right asymptotic behaviour when z or w are 
near the cusps of T. This means that we have to do some more work to find suitable bounds on 
the Green function gr r (z,w) in this case. 

Let D and D denote the open and closed unit discs in C, respectively. We recall that the 
Poisson kernel on D is defined by 



^(0 



i-ICI 2 
IWI 2 

oo oo 



n=l 



n=l 



We will use the notation 

P(t,0 = P{exp(2mt)C). 
Lemma 5.1. The Poisson kernel satishes 

P(a, ()P(-a, r])da = P((r)) for all (, n E D 



(5.1) 



and 



P(t,0 



d_ 
~dt 



t + 



2ni 



1 - exp(-27T^)C , 1 - C 



1 — exp(2irit)( 



for all (eD. 



(5.2) 



Proof. The first claim can be verified in several ways, for example using the residue theorem, 
Fourier series, or the fact that the Poisson kernel solves the Laplace equation with Dirichlet bound- 
ary conditions. The second claim is straightforward to check. □ 

Let gify denote the Green function for the Laplace operator on D; this is an integral kernel 
for the Poisson equation A/ = g with boundary condition / = on dD. It is given explicitly by 



C-v 



for all C, V <E D with C 7^ V- 



For all £ € D and t S R, we write 



1 



A(£,t) = ^— (log(l - exp(-27rztK) - log(l - exp(27rii)0). 



Lemma 5.2. The function \(£,t) satisfies 



Jo y L 



< — for all t > 0. 



Proof. We expand A(£, i) for £ e D in a Fourier series: 



£ n exp(27rmi) 



2m 

OO 

1 E 



E£" exp(-27rmi) \ 
n J 



71 = 1 



£™ sin(27rni) 
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This implies 



Jo ^ " ■ 



y 



p27rnt 



*h n Jo 



SIM 



dx 



i 00 tn 

= (si(0)-si(27rni)). 

n=l 

Here si(y) is the sine integral function normalised such that lim^oo si(y) = 0: 

si(y) = / ^^dx. 

Jy X 



It is known that 



From this we get 



7T 1 

si(0) = — and |si(x)| < — for all x > 0. 



t M^y) d i^ g"si(27r^) 



V 



2 ' n 7r 

n— 1 n=l 



1 



,iog(i-o--E 



1 ^ si(27rnt) 



n=l 



n 



and 



E 

n=l 



£" si(27rrrf) 



~ E n (2imt) 

n=l x ' 

1 °° 1 
= — V — 

2?rf ^ n 2 

n—1 

~ 2nt 6 ' 

This proves the claim. 

For S > 1 and u > 1, we write 

Js(u) = max{0, L(u) — L(S)}. 

For ^ G D, 5 > 1 and e > 0, we write 

(<t) 2 



^(0-/^1+ 2 



P(exp(27rit)C)rfi. 



Lemma 5.3. The function N$, e satisfies 



N s , e (0 arctan \/ ^r- + ir lo s I 1 - CI 

£ 7T V 2 Z7T 



< era for att £ e £>, 



where 



t\$ = 



24tt V V 5 - 1 



2 M- 1 

+ arctan 



Proof. We note that 



where 



* < T, 



V2S-2 



□ 



(5.3) 
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By the definition of J$, this gives 

NsAO = f + t^p) - L{5))P{t,Qdt. 

Using (5.2), integrating by parts, and taking the contributions for positive and negative t together, 
we obtain 



e 2 tL\l + (etf/2)(2t + A(£,t) + \(£,t))dt. 
Using the definition (1.1) of L and rearranging gives 

NsAO = I [ A - ^) (2t + A«, t) + Ml t))dt 



2^ Jo 



2(e<)2 rfi-Wy&tnKi,*)* (5 . 4) 



A+(et) 2 J 2irJ \ 4+(et) 2 J t 



i r 8 .if, ^ 



27r7 4+(et) 2 2tt 7 4 + (et) 2 t 
We consider the two integrals in the last expression one by one. As for the first integral, we have 

per/2 



,T g 2 rer/z 2 

/ -; ; — ^7"^ = - I ~dx 

Jo 4+(ei) 2 ej 1 + x 2 



4 er /r r s 

= - arctan — (5.5) 
e 2 



4 M-l 
= - arctan y — - — . 

As for the second integral in (5.4), let us write for convenience 

4 A(£,i) + A(£t) 



f 

Jo 



4 + (et) 2 t 

and 



dt 



Jo y 



Then we have 

A/ (t) = A «,«) + A(€,*) and A£(0)=log | 1 + C |. 
Integration by parts gives 

* = - lo § I 1 - « I + ITW A?(T) + I ¥TWr A ^ t)dt 

By Lemma 5.2, it follows that 

i . . I 4 1 f T 8e 2 t 1 , 

^ +1 ° g|1 - ?l ^ + 4T(^&7 + i (4+(et)»)»« dt - 
The integral can be evaluated by elementary means, and the result is 

\h + log |1 - C|| < ^ + ^ arctan y 



e / / 2 /(5-1 
+ arctan 



12 VV 5-1 

Combining this with (5.4) and (5.5) proves the claim. □ 
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Lemma 5.4. Let T be a cofinite Fuchsian group, and let S > 1 and e' > e > be real numbers 
satisfying the inequalities 

(5 + VS 2 - l) 1/2 e' < min C c ( 7 ) and (6 + ^fp^l)^ < e' . 

(a) For all z, w G H with y c {z) > 1/e' and > 1/e' and ah 7 € T, we have 

w(z, 7 u>)<<5 7 G T c . 

(b) For all z, w G H such that y c ( 2 ) > V e and such that the image of w in T\H h'es outside D c (e'), 
and for all 7 G T, we have ^(2;, 7 u>) > S. 

Proof. Let z, w and 7 be as in (a). We write 

-1 / fa b 

Suppose 7 ^ T c . Then our assumptions imply 

\c\/e'> (S+VS 2 -l) 1/2 . (5.6) 

We have 

u(z, jw) = u(a^ 1 z, a~ 1 -fw) 
= u{a~ 1 z, r y'a~ 1 w) 
= |cr- 1 Z-7 , cr- 1 w| 2 



> 1 + 



2(3?crr 1 z)(3?7'o-r 1 w) 
2(3cr c - 1 z)(97'cr c - 1 w) 

1 / 3(T^ 1 Z ^'(T^ 1 !!; 



2 \3 ; 7 / (T c w 3cr c 2; 
1 yt(z)\cv7 lw + d \ 2 + Vc( w ) 



From (5.6), we deduce 



Vc(w) yAz)\ccr c 1 w + d\< 



y^z^co-^w + d\ 2 > y c (z)(cy c (w)) 2 
Vc(w) y c (w) 

= c 2 y c (z)y c (w) 

> (MA') 2 



> 6+ VS 2 - 1. 

Using the fact that the function x ^ x + x^ 1 is increasing for x > 1, we obtain 



.<^)>i((*+v/^) + JT ^ f =) 



= 5. 

This proves (a). 

Now let z, w and 7 be as in (b). Our assumption that the image of w in T\H lies outside D^(e') 
implies 

y c ( 7 w;) < 1/e' 

and hence 

fr(g) 

y,(jw) ~ e 



>5+^5 2 -l. 
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Using the fact that the function x ^ x 1 is increasing for x > 1 as in the proof of (a), we get 

u(z, jw) = u(a~ 1 z, a^ 1 "/w) 



2 \y c {iw) y c (z) J 

>l( (J+v ^TT) + _i f =) 



= 5. 



This proves (b). 



□ 



In the following proposition, we extend our bounds on gr r to the neighbourhoods D c (e c ) of 
the cusps. We make the following abuse of notation: for z e H and S a subset of T\H, we write 
z £ S if the image of z in T\H lies in S. 

Proposition 5.5. Let T be a cohnite Fuchsian group, and let S be a reai number with 6 > 1. For 
every cusp c of T, let e' c > e c > be real numbers satisfying the inequalities 



e' c (S + VS 2 - 1) 1 < min C c ( 7 ), and (5 + \J b 2 - l)e c < e' c 
and smali enough such that the discs D c (e' c ) are pairwise disjoint. Let 

y = (t\h.)\\_\d c (c c ). 

c 

Let A and B be real numbers satisfying 

A < gr r (z, w)+ { L ( u ( z > w )) - L ( S )) < B for a11 z,w eY. 

u(z 1 *yw)<8 

(a) If c is a cusp such that z £ D c (e c ), w G F and w £ D c (e' c ), we have 

A < gr r (z, w) \- log(e c y c (z)) < B. 

voip 

(a!) If c is a cusp such that w G D c (e c ), z eF and z D c (e' c ), we have 

A < gr r (z,w) log(e c y c (w)) < B. 

voir 

(b) If c, t> are two distinct cusps such that z e D c (e c ) and w e D 5 (e ), we have 

-4 < gr r (z, w) \- log(e c y c (z)) \- log(e c y c (w)) < B. 

voir voir 

(c) If c is a cusp such that z, w € D c (e' c ), we have 

A t < g r r (z, w) - #(r n {±i}) • i- log | ge («) - 9c H| 



where ^4 C and -B c are defined using the function r$ from (5.3) by 

i c = A + #(rn{±i}) 
B t = B + #(rn{±i}) 



l 



l 



1 arctan 

7T 



1 arctan 

7T 



(5-1 



5-1 



+ e'r s 



(5.7) 
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Proof. In view of §2.1 (or Lemma 5.4(a)), the discs D c (e' c ) D D c (e c ) are well defined. Further- 
more, the assumption that the discs D c (e' c ) are pairwise disjoint implies that for every cusp c, the 
boundaries of D c (e c ) and D c (e' c ) are contained in Y. 

Let us prove part (a). We keep w € Y fixed and consider gr r (z, w) as a function of z e D c (e c ). 
The defining properties of gr r imply 

gr r (z,w) = ^j-log(e c y c (z)) + h w (z) for all z G D c (e c ), 

with h w a real- valued harmonic function on D c (e c ). By construction, h w (z) coincides with gr r (z, w) 
for z on the boundary of D c (e c ). This implies 

M 2 ) = / gr r {< T c{a + i/e c ),w)P(-a,q c (z)exp(2ir/e c ))da. 
Jo 

By Lemma 5.4(b), there are no 7 <G T such that u(z 7 jw) < S. By the assumption (5.7) on A 
and B, we conclude 

A<h w (z)<B for all z e £> c (e). 

This proves (a). Part (a') is equivalent to (a) by symmetry, and (b) is proved in a similar way. 
It remains to prove part (c). We identify D c (e' c ) with the closed unit disc D via the map 

D c (e' c ) D 

z\ — > Cz = <7c(z)cxp(27r/e' c ). 

Let grg ( ( e /) be the Green function for the Laplace operator on D c (e' c ), given in terms of gr D by 

gr 5c(ei) («,t«) = #(rn{±i})- g r D (c„C»)- 

The factor #(F n {±1}) arises because of how we defined integration on F\H in Section 1. 
Fixing w and considering gr r as a function of z, we have 

gT T (z,w) = gr Bc(e , c) (z,w) + -i-log(e' c y c (z)) + h w {z) for all z G D c (e' c ), 

where h w is a real- valued harmonic function on D c (e' c ). By construction, h w (z) coincides with 
gr r (z,u>) for z on the boundary of D c (e' c ). This implies 



h w (z)= / gr r (o- c (a + i/e' c ),w)P(-a,( z )da. 
Jo 



Applying the same argument to gr r (a c (a + i/e' c ),w) as a function of w, we obtain 

gr r (z,w) = gr Dc{<) (z,w) + ^-log(e' c y c (z)) + log(e' c y c (w)) +K(z,w), (5.8) 

where K is the function on D c (e[) x D c (e' c ) defined by 

K(z,w)= [ f gT r (o- c (a + i/e' c ),o- c (b + i/e' c ))P(-b,( w )P(-a,<; z )dbda. 

Ja=Q Jb=0 

In (5.7), we may replace T by T c in view of Lemma 5.4(a), i.e., we have 

A<gr r (z,w)+ (L(u(z,-/w)) - L(8)) < B for all z, w e dD c (e' c ). 

u(z,~/w)<5 
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Substituting this in the definition of K(z,w) and "unfolding" the action of T c , we get 

A < K(z, w) + #(r n {±1})M(C„ Cw) < B, (5.9) 
where M is the function on D x D defined by 

M(C,v)= [ f J S (u(<T c (a + i/e' c ),CT c (b + i/e' c )))P(-b,r,)P(-a,Odbda. 

Ja=0 JbeR 

Making the change of variables 

b = a + t 

and noting that 

u(a c (a + i/e' c ),a c (b + i/e' c )) = u(a + i/e' c ,b + i/e' c ) 

{b-af 



1+ 2/6? 



we obtain 



M(C,v)= [ [ J S (l+ ( -^-)p(-a-t,r ] )P(-aX)dtda. 



Interchanging the order of integration, noting that 

P{— a — t,rj) = P(a, cxp(2TTit)fj) 
and using (5.1), we simplify this to 



M{Cl V) = / r JS (* + ^2~) P[tl ^ 



Applying Lemma 5.3, we conclude from (5.9) that 



K(z, w) < B + #(r n {±1}) (J^ log 1 1 - CzCw | - j ■ \ arctan + <r s 



s 






2 


5 




1 


2 



K (z, w) > A + #(r n {±1}) (±- log |1 - C»Ct»l - 7 • | arctan y — - - ,' c r. 
We now note that 

Cz-Cu 



gr SeK) («,«;) = #(rn{±i}).i- log 



2tt 

= #(rn{±i}) ' 



i - 



^log| 9c (z) - g c (w)| + - ^ log 1 1 -Cz(u 



Combining this with (5.8) and the above bounds on K(z, w) yields the proposition. □ 
Appendix: Bounds on Legendre functions 

In this appendix, we prove a number of bounds on the Legendre functions Pj?(z) and Q£(z) that 
are used in the rest of the paper. 
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Lemma A.l. Let u € [1,3] and A > be such that X(u — 1) < |, and let s £ C be such that 
s(l — s) = A. Tiien tie reai number P s l 1 (u) satisfies the inequalities 



(2- 4/^1^ <P~_\(u) 



u- 1 

U + l' 



Proof . We start by expressing the Legendre function in terms of GauB's hypergeometric func- 
tion F(a, b;c;z). Because of the many symmetries satisfied by the hypergeometric function (see 
Erdelyi et al. [4, Chapter II]), there are lots of ways to do this. Using [4, § 3.2, equation 3] gives 



u- 1 

u + r 



F[ s,l-s;2; 



1 — u 



Next we use the hypergeometric series for F(a, b; c; z) with z < 1 (see [4, § 2.1, equation 2]): 



F{a,b;c;z) = ^ 



n=0 



(a)n(fr)n 

(c)„n! 



(A.l) 



where 



(y) n = T(y + n)/T(y) = y(y + 1) • • • (y + n - 1). 



Putting a: = for a moment, using (A.l) and applying the triangle inequality, we get the bound 



|F(s,l-*;2;-aO-l| < £ 



(s)„(l - s) r 



(2) n «! 



The assumption A(u — 1) < ^ is equivalent to Ax < |. Therefore the n-th term in the series on 
the right-hand side can be bounded as follows: 



This implies that 



(s)„{l - s) n 



(2)„n! 



< 



IlkZo(s(l-s)x + k{k + l)x) 
(2)„n! 

nro(i+fc(fc+i)) 



(2)„n! 



(2)71(2)71 



(2)„nl 



\F(s,l-s;2;-x)-l\<F(±,±;2;l)-l 
= 4/tt - 1, 

where the last equality follows from the formula 

F(a, b- c; 1) = ffife^ - a ). for 3fc> and 3fc > R(o + 6) (A.2) 
1 (c — a)I (c — 6) 

(see Erdelyi et al. [4, § 2.1.3, equation 14] or Iwaniec [10, equation B.20]) and the fact that T(3/2) = 
y/n/2. We conclude that 



P- 1 - 



u-l 



u+ 1 



u-l 



u + 



-|F(«,l-a;2;-a:)-l| 



< 



u-l 
u+ 1 



(4/7T-1) 



(4/TT-i; 



It- 1 
U+l' 



which is equivalent to the inequalities in the statement of the lemma. 



□ 
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Lemma A. 2. For all real numbers v > and u > 1, the real number Q' v {u) satisfies 



1 



u + 1 I u 2 — 1 



< QU«) < o. 



Proof. We express in terms of the hypergeometric function using [4, §3.6.1, equation 5, and 
§ 3.2, equation 36]: 



Ql(u) = 



2 y i t(i + u)t(2 + u) 



u+lj u 2 - 1 T(2 + 2u) 



u+1 



Since 2/{u + 1) < 1, the hypergeometric function is given by the series (A.l). The non-negativity 
of all the arguments gives the bounds 



0< F [v,l + v-2 + 2v; 



u + 1 

<F(v,l + v;2 + 2v;l) 

r(2 + 2i/)r(l) 
~ r(2 + v)T{l + v) ' 

the last equality follows from (A. 2). Combining this with the above formula for Q' v (u) yields the 
claim. □ 

Lemma A. 3. For all real numbers a, b and y with b > a > 0, we have 



exp 



1/1 1\\ (a 2 +y 2 ) a / 2 - 1 /4 



12 V « bj J {b 2 + y 2 ) b /2-V4 



< 



r(a + iy) 



T(b + iy) 



< exp [ b — a + 



1/1 1\\ (a 2 +y 2 ) 



,2\o/2-l/4 



12 Va &// (fe2 +y 2)6/2-l/4 



Proof. We use Binet's formula for logT (see Erdelyi et al. [4, § 1.9, equation 4]): 

/>oo 

logT(z) = (z - |) logz - z + ilog27r+ / B(t)exp(-zt)dt for 3ftz > 0, 

Jo 



where 



B(t) = 



(exp(t) - l)" 1 - t- 1 + \ 



| coth | - 1 



f 2 



We write 



/>oo 

M(a, y) = $l B(t) exp(-(a + 
Jo 

exp(— ai) cos(yt)dt. 



Then we have 



log 



r(a + zy) 



r(6 + iy) 



((a + - 1/2) log(a + iy) -(b + iy- 1/2) log(6 + iy)} 



We note that 



- a + b + M(a, y) - M (b, y) 

= (a - 1/2) log \a + iy\ - (b - 1/2) log \b + iy\ 

- y arg(a + iy) + y arg(6 + iy) - a + b + M(a, y) - M(b, y). 



(b — a)y 

y arg(a + iy) - y arg(6 + iy) = y arctan — - € [0, - a\. 



ab + y 2 
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Using b > a, we conclude 
T(a + iy) 



log 



and 



T(b + iy) 

r(o + %) 



< (o - 1/2) log \a + iy\ - {b - 1/2) log \b + iy\ - a + b + M(a, y) - M{b, y) 



log 



> (a - 1/2) log |a + iy| - {b - 1/2) log |6 + iy\ + M(a, y) - M{b, y). 



T{b + iy) 

It remains to bound M(a, y) — M(b, y). The function B satisfies 
< Bit) < lim Bix) = 1/12 for all t > 0. 

Using this and the positivity of exp(— at) — cxp(— fei), we bound M(a, y) — M(b, y) as follows: 

/>OC 

\M(a,y) - M(b,y)\ < / B(i)(exp(-ai) - exp(-bt)) \cos(yt)\ dt 
Jo 

= B(t)(exp(-at) - exp(-bt))dt 
Jo 

1 f°° 

< — (exp(— at) — cxp(—bt))dt 
12 Jo 

12 \a b) 

This implies the inequality we wanted to prove. 

Corollary A. 4. For b > a > 0, b > 1/2, and y <G R, we have 



□ 



r(a + iy) 



T{b + iy) 



< exp [b — a 



1 /l 1 



12 V a b 



(a z + y z ) 



2\-(6-a)/2 



Given a real number cr e (0, 1/2), we consider the strip 

S„ = {s e C | <7 < 5Rs < 1 - a}. 

We put 



C a = max{l,tan7T0-}((T 1 — 1) 1//4 exp ( - 



2 24a(i+a) 



= max{l, tan7rcr}(cr 1 — 1) 1//4 exp [ i 



2 24(l-a)(|-<7) 



(A.3) 



Proposition A. 5. Let m be an even non-negative integer, and let a G (0, 1/2). For all s e S a 
and all u > 1, we have 



|p-_7H|<| s (i- s )| 



-(2^+1)74^2" + C' a x m - 1+ ° ^ |(|-m) n | _ 2n 



0F(4(m 2 - 1))W2 



E 

n=0 



where 



x = u + \]u 2 — \, 



x + x 



Proof. We use the following expression for P s ™[ (see Erdelyi et al. [4, §3.2, equation 27]): 



P s -- m Au) 



r(-i + s ) x 



m+s— 1 



v / 7rr(m + 5) (x - x- 

r(i- s ) 



— -F(±-m,l-m-s;§-s;a: 2 ) 



v / 7rr(m + 1 - s) (x- x- 1 ) 
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l\m r V2 



m, —m + s; 7} + s; x 2 ) . 



Using the hypergeometric series (A.l) and the functional equation 

T(z)r{i-z) = ^—, 

smwz 

we get 

T(m + s) ^ (f-*)n n\ 



- S ) m-s y (I ~ m )»(~ m + f)n X 2 " 

+ r(m+l-,s) X ^ (I + S )„ n! 



OO / l 



(i- m ) n j-r(-l + S ) T(l- S )T(n-m+l- S ) m _ 1+s _ 2n 



^ n! t r ( m + s ) r(n + | - s)r(-m + 1 - s) 



r(i-.s) r(± + s )r(n~ m + s ) _ 2n ' 

r(m + 1 - s) T(n + \ + s)T(-m + s) 



(5 — TO )n f sin 7r(m + s) T(n — m + 1 — s) 



n=0 



,m— 1+s— 2n 



sin7r(-i + s) T{n+\-s) 

sin7r(m + 1 - s) T(n - m + s) ^ m _ s _ 2 »| 



sin7r( 2 - — s) r(n + | + s) 
Basic trigonometric manipulations simplify this to 



(-l) m tan(7T S ) ^ (£-m) n f I> - m + S ) 
- lW Vi(x-x-i)-^ n! 1 



jn—s—2n 



V^(x - x- 1 )™ ^ n\ lr(n+i + s) 

(A 4) 



T(n-m+l-s) ^ m _ 1+s _ 2ll 
r(n+f -*) 



-:r" 



On the right-hand side, the pole of tan(7rs) at s = 1/2 is cancelled by a corresponding zero of the 
function defined by the sum. 

For fixed u > 1, we consider the holomorphic function 

irOO = (*(i - .)) (2m+1)/4 tanM £ ii-^i { + ro _ s _ 2 „ 

_ r(n-m + l-s) ^ m _ 1+s _ 2ra ] 
3 



T(n+l-s) 



on where we have fixed a branch ofm (s(l - s))< 2ro+1 )/ 4 . Because ff(s) = - s), the 
Phragmen-Lindelof principle gives 

sup |-ff(s)| < sup \H{a + iy)\. 
ses„ yen 

Together with (A. 4), this implies 

|s(l - s)r (2m+1)/4 su Pt , eR \H(a + iy)\ 

Let y G R and s = a + iy. Then we have 

- s )|( 2m+1 )/ 4 = [a 2 +? /)(2™+l)/8 ((1 _ ct) 2 +y 2 )(2m +l)/8_ 

A straightforward calculation gives 

|tan7rs| = |tan7r(<7 + iy)\ < max{l,tan7r<7}. 
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Using Corollary A. 3 and the assumption that m is even, we bound the quotients of T-functions 
appearing on the right-hand side of (A. 5) independently of n: 



Tin — m + s) 



I>+± + ,s) 



E>+£)j 
|r(n+| + a)| . 11 \j + s\ 



n — 1 1 

n u4- 



<exp(--\ J- = -)\n + a + iy\ 1/2 — 



+ iy\ m / 2 \l - a + iy\ m / 2 



< exp - + 



2 24ct(± + a) J {a 2 + y 2 )(™+i)/4((i - a ) 2 + y 2 ) m / A ' 



This implies 



\s{l-s) 



,(2m+l)/4 



r(n — m + s) 



r(n + i + s ) 



< exp - + 



2 24a(i+a) 



: exp 



(a 2 + ^)(2r»+l)/8( (1 _ g )2 + jj 2 ) (2m+l)/8 
((J 2 + j/2)(m+l)/4(( 1 _ (T )2 + y 2yn/4 
1 1 \ ((1-^)2+^)1/8 



< exp - + 



2 ' 2Aa{\+a)) (a 2 +y 2 )i/8 
1 1 \ (1-a) 1 / 4 



and hence 

Similarly, 

r(n — m + 1 — s) 



I , .i(2m+l)/4, 

s(l - s) |tan7rs 



2 2Aa{\+a)J cr 1 ^ 
r(n — m + s) 



r(n + i + s ) 



< a. 



T(n+l-s) 



< 



exp 



2 24(1 -a){%-a)J ((1 - ct) 2 + y 2 )(™+i)/4( CT 2 + y2 ) m ;" 



and 



9(1 - S) 



I (2m+l)/4 



r(n — m + 1 — s) 



T{n+l-s) 



< exp - 



2 24(1 - <7)(§ - a) 



(a 2 + y 2 )( 2m + 1 )/ 8 ((l - a) 2 + y 2 )(2m+i)/8 

((1 - a) 2 + y2)m/4( CT 2 + y 2)(m+l)/4 
1 1 \ ((1-^+^)1/8 



This implies 



We conclude that 



< exp 



| s(1 _ s)|(2m+ l)/4 |tan7rs| 



" X1)( 2 + 24(1 -*)(§- a),/ (a 2 +y2)i/8 
1 1 \ (1-ct) 1/4 



2 24(l-<7)(|-a); a 



rV4 



r(n — m + 1 — s) 



r(n+f -s) 



sup |ff (a + »i/)| < £ K ' ^^ ) " l (C ff a;" t - g - 2 " + C> m - 1+CT - 2 ") 



sea 



n=0 



n! 



00 i c l \ I 



n=0 



n! 



-x 



This hnishes the proof. 



□ 
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With C a and C' a as in (A. 3), we define an elementary function p a : [l,oo) — > R by 



Pg ( U ) = ^ " +C'° xl+a ((1 _ ^-2)3/2 + 3^-2-) ) where ^ = u + ~ x (A 6) 

Corollary A. 6. For all a e (0, 1/2), s e S a and u > 1, we have 



i*7- 2 i(«)i<Ki-*)r 6/4 



m 2 - 1 

Proof. We note that 

l(-|)n| 



2jn| •(-§)„ otherwise. 



This implies that for z € (0, 1), we have 

E K~f )"l n _ (~|)" n _ g (~|) 1 1 

n! ^ n! 1! 

n=0 n=0 

= (l-z) 3 / 2 + 3z. 

The claim immediately follows from this identity and Proposition A. 5. □ 
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